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The two-orbital Hubbard model with the Hund coupling is investigated in a metallic phase 
close to the Mott insulator. We calculate the one-particle spectral function and the optical 
conductivity within dynamical mean field theory, for which the effective impurity problem is 
solved by using the non-crossing approximation. For a metallic system close to quarter filling, a 
heavy quasi-particle band is formed by the Hubbard interaction, the effective mass of which is 
not so sensitive to the orbital splitting and the Hund coupling. In contrast, a heavy quasi-particle 
band near half filling disappears in the presence of the orbital splitting, but is induced again by 
the introduction of the Hund coupling, resulting in a different type of heavy quasi-particles. 
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§1. Introduction 

Strongly correlated electron systems, such as the 
high-Tc -superconductorall-' and the heavy fermion com- 
pounds p-* have attracted much interest. In order to deal 
with elefjtron correlations, the single-orbital Hubbard 
modelBQEP has been studied extensively with various 
analytical and numerical methods. ^-Among others, dy- 
namical mean field theory (DMFT)® enables us to sys- 
tematically investigate strong correlation effects. By us- 
ing DMFT characteristic properties related to the Mott 
transition has been described successfully for the single 
orbital casccP l-. 

In real materials such as the transition metal oxides,tP 
orbital degeneracy plays an important role, giving rise to 
a wide variety of interesting phenomena. In such cases 
the Hund coupling and the orbital splitting are the addi- 
tional key parameters to discuss electron correlations. It 
is thus desirable to study the orbitally degenerate Hub- 
bard model in order to better understand strongly corre- 
lated electron systems. Recently, a number of theoretical 
studies on the orbitally degenerate Hubbard model have 
been dope .by. using the slave-boson mean field approx- 
imation jlJ'EilJl3.^Jic_)^^ Monte Carlo methodtl^ 

and DMFT.yyy.By' 

In this paper, we systematically study electron correla- 
tions for the two-orbital Hubbard model with particular 
emphasis on the interplay of the Hubbard interaction, 
the Hund coupling and the orbital splitting. By em- 
ploying DMFT combined with the non-crossing approx- 
imation (NCA), we compute the one-particle spectral 
function and the optical conductivity. We discuss how 
the heavy quasi-particle behavior shows up in a metallic 
state close to the Mott insulator with orbital degeneracy. 
We find that the origin of heavy quasi-particles can be 
different between quarter fiUing and half filling. In par- 
ticular, we point out that the Hund coupling plays an 
essential role to form a heavy quasi-particle band in the 
vicinity of half filling. 

This paper is organized as follows. In the next section. 



we briefly describe the model and the method, and then 
in §3 we show the results obtained for the one-particle 
spectral function and the optical conductivity. A brief 
summary is given in §4. 

§2. Model and Method 

2.1 Two- Orbital Hubbard Hamiltonian 

We study the Hubbard model with two orbitals. The 
Hamiltonian is given by 

"2" 



H = 



E 



E 



.mm' t 



U' ^ nimT"-imi + E '^iicr'^i^S- 
i.m 2,(7 
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J 



(2.1) 



where CimaiCma) annihilation (creation) operator 

of electrons with the orbital m{— 1, 2) and spin fT(=t, i) 
at the site i and A represents the energy splitting be- 
tween the two orbitals. Here, U'{U) is the Coulomb re- 
pulsion between electrons in the same (different) orbital 
and J is the ferromagnetic Hund coupling. When the 
two orbital-levels are degenerate (A = 0), the system is 
rotationally invariant with respect to the spin and orbital 
degrees of freedom, leading to the condition, [/' = [/ -j- J. 
We use this relation and assume that the hopping is diag- 
onal with respect to the orbital indices (i™™ = —tSmm')- 
We shall deal with a paramagnetic metallic phase close 
to the Mott insulator in this paper. 

2.2 Dynamical Mean Field Theory 

It is known thatllMFT is justified in the limit of large 
spatial dimensionsErLP and gives a rather good approxi- 
mation even in three dimensions. In this method, a lat- 
tice electron problem is replaced by the corresponding 
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impurity one embedded in an effective medium deter- 
mined seif-consistently. Tfie effective local Hamiltonian 
is given by 



Hcff — H\oc + i?mcd, 



(2.2) 



We now introduce the ionic resolvent and its self- 
energy, which are represented in terms of the basis set 
I a) as, 

R^[u)=[uj-E^-^^{lo)]-\ (2.13) 



U ^ ni^n2a + ([/ - J) E '^1^ 



n2a 



J{c\^ciic\^C2-^ + h.c), 



(2.3) 



h.c.) 



E 



(2.4) 



where the impurity level in the reduced problem is 



E 



fm 



Ef + /\/2 (orbitall) 
Ef~/\/2 (orbital 2) 



(2.5) 



Here, Vkma represents the hybridization between the im- 
purity site and an effective medium. 

Let us first introduce the so-called cavity Green func- 
tion and the full Green function, which are defined re- 
spectively as 

Qo,na[uj) = [CJ + i(5 + /i - S/ni - A™<^(w)]"^ (2.6) 

Gma{u:) = [QoL{^) - ^ma{i0)]-^ (2.7) 

where /i is the chemical potential. Here A,„cr(a-') repre- 
sents the hybridization function given by 

In DMFT, the self-consistency equation is given by 

[lo + i5 + ii- Ef,n - ^ma{uj) - — G^^ (w)] -\(2.9) 

where we have employed the semi-elliptic density 
of states for non-interacting electrons, iVo(e) = 
2/{TrD'^)VD^ - £2 with the bandwidth D. 



In order to solve 
further use NCaM 



t jvn im purity problem, we 



feSSk^) which is known to 
provide rather reliable results ip, agreement with the 
quantum Monte Carlo calculation^ for the single-orbital 
case. We first diagonalize the local Hamiltonian iJioc, 



16 



H\oc = E Ea\a){a\, 



(2.10) 



where \a) and Ea denote the spin-orbital eigenstate and 
the corresponding energy for the isolated ionic system. 
Using these eigenstates, the fermionic annihilation oper- 
ator is represented as 



E^-'i")("'i' 



(2.11) 
(2.12) 



X / def{T]aa'e)Amaie)Ra'ii^ + Vaa'e), (2.14) 
J oo 

with the Fermi distribution function /(e), where r]aoi' 
takes 1(— 1) if the particle number of the state \a) is 
larger (smaller) than that of \a'). In eq.(2.14), we have 
neglected the vertex corrections, as often done in the 
treatment of DMFT. The full Green function is now 
given by 

G,nA^) = J2 (2.15) 

aa' 

where Gaa' (i^) is written in terms of the ionic resol- 
vents, 

Go„'(c^)-| I dee-'3^ 
X [paie)Ra,{Lu + e) ~ p^,{e)Ra{e - uj)] (2.16) 

and 

^ = E /dee-%a(e). (2.17) 

In eq. (2.16) and (2.17), pa(w) = -l/7rlm(i?a(w)) is the 
spectral function of the ionic resolvent. 

This completes the self-consistent procedure for 
DMFT with NCA. 

§3. Numerical Results 

We numerically iterate the procedure in the previous 
section until the calculated quantities converge within 
desired accuracy. In the following discussions, the band- 
width D is taken to be unity for simplicity. 

3. 1 Heavy quasi-particles near quarter filling 

Let us start with a metallic system close to quarter 
filling (ritot ~ 1)- We calculate the one-particle spectral 
function defined by 



Pn 



r{uj) = ImG,, 



(3.1) 



for several different values of the Hubbard interaction U. 
We first discuss the case in the absence of the Hund cou- 
pling, J — 0. The results are shown in Fig.l. As the Hub- 
bard interaction increases, a sharp peak structure ap- 
pears and gets narrower around the Fermi level, 
caused by the formation of heavy quasi-particles Jlil^E^ I 
At the same time, it is seen that a pseudo-gap structure 
develops just above the Fermi level because the Hub- 
bard interaction suppresses the charge fluctuation while 
enhances the spin fluctuation near quarter filling. This 
behavior is similar to that expected for the hole-doped 
Hubbard model with single orbital. The other peaks are 
easily identified with the singly, doubly and triply oc- 
cupied states, respectively. The intensity of the fully 
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occupied state is too small to observe in this figure. 

The formation of a heavy quasi-particle band is more 
clearly seen in the quasi-particle renormalization factor 
shown in the inset (b) of Fig.l, which is defined by 



Z 







duj 





(3.2) 



Since the mass enhancement factor of quasi-particles is 
inversely proportional to the renormalization factor, it 
is confirmed that the mass of the renormalized band in- 
deed increases as the Hubbard interaction U increases. 
In Fig. 2, we show how the Hund coupling J affects 
the spectral properties near quarter filling. Although 
the spectral weight for the doubly occupied state splits 
into two parts as J increases, the quasi-particle spectrum 
near the Fermi level is not so much affected by weak or 
intermediate J. In particular, for the strong Hubbard 
interaction, the Hund coupling hardly affects the quasi- 
particle mass although it modifies the spectrum in the 
high energy regime. 

Similarly, even when the splitting between two orbitals 
is introduced, the formation of quasi-particles is not so 
much affected, as seen from Fig. 3. For the lower 
band, the spectral weight is gradually transferred from 
the doubly and triply occupied states to the singly occu- 
pied state as A increases, while for the upper band, the 
weight for the singly occupied state decreases, as should 
be expected. Note that a quasi-particle band always ex- 
ists in the lower band irrespective of A. We can check 
that the renormalization factor for the lower band is al- 
most unchanged as a function of A. (Fig. 3, inset) For 
sufficiently large A, the spectrum of the lower band is re- 
duced to that of the single-orbital model near half filling, 
which is indeed seen in the case of A = 0.8. 

Summarizing, the formation of heavy quasi-particles 
near quarter filling is caused by the Hubbard interaction, 
as is the case for the single-orbital model. In particular, 
a heavy quasi-particle band and a pseudo-gap structure 
are almost always formed for reasonably large U, irre- 
spective of the strength of the orbital splitting and the 
Hund coupling. As we will see momentarily, this is not 
the case for the formation of heavy quasi-particles near 
half filling. 

3. 2 Heavy quasi-particles near half filling 

We now turn to a metallic system close to half-filling 
(ntot ~ 2). Let us first discuss the system without or- 
bital splitting. Fig. 4 shows the one-particle spectrum as 
a function of the Hund coupling J. As mentioned above, 
a heavy quasi-particle band exists around the Fermi level 
for [7 = 2. In contrast to the quarter filling case, the 
spectrum near the Fermi level is considerably affected by 
the Hund coupling, since the energy eigenvalues of two- 
particle configuration are quite sensitive to the strength 
of the Hund coupling. With the increase of J, the heavy 
quasi-particle band disappears and the charge excitation 
gap develops around the Fermi level. This implies that 
the Hund coupling has a tendency to drive the system 
to the Mott insulator together with the Hubbard inter- 
action,E3 in contrast to the quarter-filling case. 



We next discuss how the orbital splitting A affects the 
formation of heavy quasi-particles near half filling. We 
present the results for the spectral function near half- 
filling in Fig. 5. The overall structure in the spectrum has 
the shape expected naively. Namely, with the increase 
of A, the weight of the middle peak (doubly occupied 
state) gets larger and that of the top peak (triply occu- 
pied state) becomes smaller in the lower band case (a), 
while in the upper band case (b) , the opposite behavior is 
observed. It should be noted here that the quasi-particle 
band disappears gradually with the increase of A and is 
absorbed into the doubly occupied state. This tendency 
is different from that observed for quarter filling. Near 
half-filling, both of the spin fluctuation and the charge 
fluctuation are suppressed in the presence of A, since the 
lower band is almost fully occupied by electrons of spin 
singlet. As a result, the system behaves like a renor- 
malized band insulator with both of the spin and charge 
gaps. Therefore, the formation of heavy electrons is sup- 
pressed by A in a metallic phase close to half filling. 

The Hund coupling, however, gives rise to a dramatic 
change in the spectrum in this case. Fig. 6 shows the 
spectral function for A — 0.8 in the presence of the Hund 
coupling. Although a quasi-particle band around the 
Fermi level is almost invisible in the absence of the Hund 
coupling (J = 0), it is enhanced again with the increase 
of J in both (a) and (b). In the absence of the Hund cou- 
pling, each lattice site is mostly occupied by two electrons 
of spin singlet in the same orbital. When the Hund cou- 
pling is introduced, the energy of the spin-triplet state 
is lowered, and can be degenerate with the singlet state. 
This induces the large spin fiuctuation, enhancing the 
quasi-particle mass to form a heavy quasi-particle band. 
This tendency is clearly seen in the renormalization fac- 
tor shown in Figs. 5 and 6. As seen in the inset of Fig. 
5, when the orbital splitting increases, the renormaliza- 
tion factor for both of the lower (open triangle) and up- 
per bands (filled triangle) increases, while in the inset 
of Fig. 6 the introduction of the Hund coupling again 
reduces the renormalization factor. These results imply 
that the spin fluctuation caused by the Hund coupling 
enhances the quasi-particle mass in a metallic state near 
half filling, which is not observed near quarter filling. 

3.3 Optical Conductivity 

In order to further investigate the dynamical response, 
we calculate the optical conductivity. We employ the 
following formula for the optical conductivity. 



cr(aj) = vr 



m.(T 



de' 



deiVo(e) 



X Ama[e,£)Am<T[e,uj + e) ,(3.3) 



A„icr(e,t^) = --Imf — — 



(3.4) 



where the vertex correction is neglecteil- which is jus- 
tified in the limit of large dimensions.Ll'Ej'EP The com- 
puted results are shown in Fig. 7. In contrast to the 
one-particle spectrum, the overall structure of the op- 
tical conductivity is not so sensitive to the parameters 
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employed. Nevertheless, we can see some distinct behav- 
iors near quarter filling and half filling. We first discuss 
the orbitally degenerate case (A — 0) near quarter fill- 
ing and half filling shown in Fig. 7 (a) and (b). It is 
clearly seen both in (a) and (b) for J = that the opti- 
cal conductivity consists of the Drude-like part (w ~ 0) 
and the interband part. With the increase of J in (a), 
the weight for the interband absorption gradually shifts 
to the lower-energy regime, reflecting the splitting of the 
doubly occupied state due to the Hund coupling, as al- 
ready shown in Fig. 2 (a) . On the other hand in (b) the 
Drude-likc part gets narrower and the interband absorp- 
tion shifts to the higher-energy part with increase of J . 
This result in (b) reflects that the Hund coupling has a 
tendency to enhance the charge excitation gapJlS* which 
is different from the case near quarter flUing. 

We next discuss the effect of the orbital splitting by 
observing Fig. 7 (c) and (d), where we assume that the 
Hund coupling is absent for simplicity. As seen from (c), 
near quarter flUing, the amplitude of the conductivity be- 
comes smaller with the increase of A although the overall 
structure is almost unchanged. For sufficiently large A, 
the amplitude becomes about half of the case of A = 
since the lower band corresponds to that of the single- 
band—model and the upper one exists above the Fermi 
levelEiP (see Fig. 3). On the other hand, the Drude- 
like peak in (d) becomes narrower while the interband- 
absorption weight gets smaller as A increases because 
both of the charge and spin fiuctuations are suppressed 
for large A near half filling, as seen in Fig. 5. Therefore, 
a slight hole-doping into half-filled band brings about the 
sharp Drude-like absorption peak in the low energy part. 

Finally, we discuss the interplay of the Hund coupling 
and the orbital splitting. As seen in Fig. 7 (e), the 
conductivity does not change its characteristic behavior 
except that the interband absorption slightly shifts to 
the higher energy regime with the increase of J. This 
is because the Hund coupling hardly affects the low- 
energy spectrum for large A near quarter filling, since 
the lower band is reduced to that of the single-band Hub- 
bard model as mentioned in Fig. 5. Near half filling 
(Fig. 7(f)), however, a rather drastic change occurs in 
the conductivity: the Drude-like part gets wider and the 
interband absorption is induced again as J increases, re- 
fiecting the appearance of a quasi-particle state near the 
Fermi level as shown in Fig. 6. 

§4. Summary and Discussions 

We have studied how the interplay of the Hubbard 
interaction, the Hund coupling and the orbital split- 
ting affects the formation of heavy quasi-particles in 
the two-orbital Hubbard model. In order to investi- 
gate dynamical quantities in the whole energy region, 
we have exploited dynamical mean field theory with the 
non-crossing approximation. We have discussed the one- 
particle spectral function and the optical conductivity 
for a metallic system in the vicinity of the Mott insula- 
tor. In particular, similarity and difference in the for- 
mation of heavy quasi-particles between quarter filling 
and half filling have been clarified. Near quarter filling, 
hole-doping into the Mott insulator leads to the forma- 



tion of the heavy quasi-particle band, which is slightly 
rcnormalizcd by the Hund coupling or the orbital split- 
ting. On the other hand, near half filling the orbital 
splitting prevents the formation of heavy particles, driv- 
ing to the system to a correlated band-type insulator, for 
which a heavy quasi-particle band is not formed by hole 
doping. If the Hund coupling is switched on, however, 
heavy quasi-particles appear again due to the enhanced 
spin fluctuation. 

In connection with the experimental results for transi- 
tion metal oxides such as Lai-ajSrajTiOap-* the quarter- 
filled |_madaL (ntot ~ 1) has been studied exten- 
sivelyEj^cStS' A heavy quasi-particle band studied here 
was pointed out in a different method, and the overall 
structure of the one-particle spectrum was quantitatively 
discussed in comparison with the photoemission exper- 
iments,E2P although a sharp quasi-particle band in the 
spectrum was not clearly observed experimentally. On 
the other hand, the half-filled Hubbard model with two 
orbitals has not been studied systematically thus far in 
comparison with the quarter-filled model. A possible ex- 
ample in transition metal oxides-iielevant for this case 
[ntot ~ 2) may be Lai_2;Sr2;V03EiHiP with electron con- 
figuration of 3c?^ in V ions. In this material, the Hund 
coupling forms the S ~ \ spin at each site, resulting in 
the antiferromagnetic ground state at half filling. Ac- 
cording to the results in §3.2 (Fig. 4), when holes are 
doped into this system to realize a paramagnetic jinetal, 
a heavy quasi-particle band may not be observedEir be- 
cause of a rather large Hund-coupling. Nevertheless, we 
think that for related materials where the orbital split- 
ting due to a crystalline fleld becomes comparable to the 
Hund coupling, a heavy quasi-particle formation may be 
more easily observed experimentally, reflecting the inter- 
play of the above two effects, as pointed out in Fig. 6. 
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Fig. 1. One-particle spectral function pmcr{i^) at T = 0.1 for J = 
and A = 0. The electron density per site is rttot = 0.96 (i.e. 
close to quarter filling). The energy is measured from the Fermi 
level u! = 0, and the bandwidth D is taken to be unity. The inset 
(a) shows the magnification of the low-frequency part, while (b) 
shows the quasi-particle renormalization factor Z as a function 
of U for ntot = 0.96 (open circles), ntot = 0.98 (filled circles) 
and ntot = 1-00 (triangles). 




Fig. 2. One particle spectral function for various J at T = 0.1 
and [7 = 2. The electron concentration is 0.96. The insets (a) 
and (b) show the magnification of the low-frequency behavior 
and the renormalization factor. 
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Fig. 5. One-particle spectral function for various values of A: 
U = 2, J = and T = 0.1. The system is close to half fill- 
ing, ntot = 1.90. (a) and (b) represent the spectrum for the 
lower band and the upper band, respectively. The inset shows 
the renormalization factor for lower band (open circles) and up- 
per band (filled circles). 



Fig. 4. One particle spectral function near half filling, ntot = 
1.90: C7 = 2, A = and T = 0.1. 
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Fig. 6. Onc-particlc spectral function at (7 = 2, A = 0.8 and 
T = 0.1. The electron concentration is 1.90 (near half filling), (a) 
and (b) represent the spectrum for the lower band and the upper 
band, respectively. The inset shows the renormalization factor 
for lower band (open circles) and upper band (filled circles). 
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Fig. 7. Optical conductivity o-{u!) for (a) A = and ntot = 0.96, 
(b) A = and rttot = 1.90, (c) J = and ntot = 0.96, (d) J = 
and ntot = 1-90, (e) A = 0.8 and ntot = 0.96, and (f) A = 0.8 
and ntot = 1.90. The other parameters are U = 2 and T = 0.10. 



